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Abstract 

The magnetic properties of various iron pnictides are investigated using first-principles pseu- 
dopotential calculations. We consider three different families, LaFePnO, BaFe2Pn2, and LiFePn 
with Pn=As and Sb, and find that the Fe local spin moment and the stability of the stripe-type 
antiferromagnetic phase increases from As to Sb for all of the three families, with a partial gap 
formed at the Fermi energy. In the meanwhile, the Fermi-surface nesting is found to be enhanced 
from Pn=As to Sb for LaFePnO, but not for BaFe2Pn2 and LiFePn. These results indicate that it 
is not the Fermi surface nesting but the local moment interaction that determines the stability of 
the magnetic phase in these materials, and that the partial gap is an induced feature by a specific 
magnetic order. 

PACS numbers: 71.15.Mb, 71.20.-b, 75.25. +z, 71.18.+y 



The iron pnictide superconductors and their fascinatingphysical properties have become 
central issues in many fields since their recent discoveries l|, y, |3| . The prototype materials 
are REFeAsO with variouof RE (rare earth) elements, and the superconducting transition 
temperature (T^) is as high as 55 K in doped SmFeAsO {4]. Other compounds with various 
.ypes of insulating layers are also superconducting when doped, such as K-doped BaFe2As2 



io|,lii|. 



, l6| and SrFeaAsa [7|, [Sj with of 38 K, and LiFeAs with of 16 K Q or 18 K 
Without doping, these materials exhibit a peculiar magnetic structure of a stripe-type an- 
tiferromagnetic (AFM) spin configuration coupled to orthorhombic atomic structure, and 
either hole or electron doping destroys the AFM and the superconductivity emerges subse- 
,ue„t., Heuce thejsaaetis,„,s considered to be e.ose.y related to the supereonductivty 



in these materials 



12 
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y,ii5i,ii6(, and the _spin-fiuctuation-mediated superconductivity 



is assumed in many theoretical works 
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19|. 



Understanding the nature of magnetism in these materials is thus of crucial importance. 



)ut it still under debate. On one hand, many theoretical 
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22, 
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25, 



22l | and experimental 



26l | works emphasize on the itinerant nature of the magnetism of the spin 
density wave (SDW) type, since the electron and hole Fermi surfaces (FS) are separated 
by a commensurate nesting vector in iron pnictides, which is further supported by the 



reduced magnetic moment of about 0.3 fiB 2J, |25|| and the energy gap near the Fermi 



energy {Ep). On the other hand, there are also in terp retations based on the Heisenberg-tj^e 
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291 ]. In this localized- moment picture. 



interaction between localized spin moments |27l . 
the observed stripe-type AFM ordering results from the frustrated spin configuration with 
the next-nearest-neighbor exchange interaction ( J2) larger than half of the nearest-neighbor 
(NN) interaction (Ji). The itinerant and the local-moment pictures are based on different 
assumptions on the electron itinerancy, but a more comprehensive mechanism might be 
discovered by combining the two pictures 
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Recently, motivated by the great success of the As substitution for P in LaFePO on raising 
Tc, hypothetical iron antimonide compounds have been studied as candidates for a higher-Tc 
superconductor by first-principles calculations [32,|33|. In these works, Sb substitution for As 
is found to modify the FS nesting and the magnetic stability significantly. Thus, with more 
variation of compounds including antimonides, more comprehensive understanding of the 
nature of magnetism in iron pnictides would be possible through a systematic comparative 
study dealing with many different types of compounds altogether. 
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In this study, we present our density-functional pseudopotential calculations of the 
electronic and magnetic properties of various iron arsenides and antimonides: LaFePnO, 
BaFe2Pn2, and LiFePn (Pn=As and Sb). We find that there is no systematic trend of FS 
nesting feature between arsenides and antimonides, whereas the stability and the local Fe 
spin moment of the magnetic phase increase from arsenides to antimonides for all three 
types of compounds. This finding is consistent with Heisenberg-type interaction picture 
that the local Fe moment is larger for antimonides with the enhanced Hund's rule coupling 
due to their larger lattice constants. We also find that the FS reconstruction and the sub- 
sequent formation of a partial gap in the density of states (DOS) at Ep can be regarded as 
a secondary effect caused by the magnetic ordering of local moments. 

Our first-principles calculations are based on the density-functional theory (DFT) within 
the generalized gradient approximation (GGA) for the exchange-correlation energy func- 
tional l3J] and the ab-initio norm- conserving pseudopotentials as implemented in SIESTA 
code [35]. Semicore pseudopotentials are used for Fe, La, and Ba, and electronic wave func- 
tions are expanded with localized pseudoatomic orbitals (double zeta polarization basis set), 
with the cutoff energy for re al sp ace mesh of 500 Ry. Brillouin zone integration is performed 



by Monkhorst-Pack scheme 



36| with 12 X 12 X 6 k-point grid. 



First we obtain the optimized cell parameters and atomic coordinates of compounds by 
total energy minimization, as listed in Table I. For the non-magnetic (NM) phase, tetragonal 
structures are obtained while the stripe-type AFM phase prefers the orthorhombic structure 
of the approximate x \/2 supercell, in agreement with experiments. The lowering of 
the total energy per Fe atom in the stripe-type AFM phase in the optimized orthorhombic 
structure relative to the NM phase in the optimized tetragonal unit cell is 354 and 706 meV 
for LaFeAsO and LaFeSbO, 297 and 745 meV for BaFe2As2 and BaFe2Sb2, and 153 and 
523 meV for LiFeAs and LiFeSb, respectively. Along with the local magnetic moments on 
Fe atoms displayed in Table I, this result implies the existence of a universal trend that 
the magnetism is stronger for antimonides than for arsenides irrespective of the detailed 
material properties. 

Figure 1 shows the calculated FSs on the = plane. To facilitate the investigation of 
the nesting feature, the electron and hole surfaces are drawn together in the reduced Brillouin 
zone for the x ^/2 supercell. LaFeSbO shows an enhanced nesting between the electron 
and hole surfaces which coincide with each other very isotropically with almost circular 
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shapes compared with LaFeAsO 32] • For BaFe2Pn2, the arsenide exhibits a moderate 
nesting feature, while nesting looks poor for the antimonide because hole surfaces, which 
are present in the arsenide, are missing so that the electron surfaces have no hole surfaces 
to couple with nearby. LiFeSb also shows an inefficient nesting compared with LiFeAs with 
some hole surfaces missing around the F point. 

The nesting feature can be more quantitatively estimated by evaluating the Pauli sus- 
ceptibility Xo(ct) as a function of the momentum q in the static limit with matrix elements 
ignored. The result is displayed in Fig. 2. For LaFePnO, Xo is larger for LaFeSbO for entire 
range of q, especially at the nesting vector q = (vr, tt) where the pronounced peak is located. 
This peak indicates the enhanced FS nesting for LaFeSbO, consistently with the FS topology 
in Fig. 1. For BaFe2Pn2, situation is drastically different. Although the susceptibility for 
BaFe2As2 has similar q dependence with those for LaFePnO, the susceptibility for BaFe2Sb2 
is larger only for partial range of q with very weak q dependence and moreover there is no 
peak at q = (tt, tt). This feature clearly reflects the poor FS nesting in BaFe2Sb2 due to the 
lack of hole surfaces, as shown in Fig. 1. Finally, LiFeSb also has smaller Xo(ci) than LiFeAs 
near (7r,7r), hence LiFeSb has less effective FS nesting at (tt, tt) than LiFeAs. 

Although many previous studies suggest the itinerant magnetism in iron pnictides that the 
stripe-type AFM is the SDW type driven by the FS nesting, our results are in contradiction 
with this picture of magnetism. As we have just discussed, the FS nesting for q = (7r,7r), 
at which the stripe- type AFM occurs, is more pronounced for LaFeSbO than LaFeAsO, 
while BaFe2As2 and LiFeAs have more effective nesting feature than BaFe2Sb2 and LiFeSb, 
respectively. Thus, there is no universal trend in the FS nesting feature between arsenides 
and antimonides, which is in contrast, however, with the result that magnetism is stronger 
for antimonides than the respective arsenides for all three types of iron pnictides, with larger 
energy differences between AFM and NM states and greater Fe local magnetic moments for 
antimonides. This implies that the contribution of itinerant electrons to the magnetic energy 
and moment is relatively small. 

In order to obtain a deeper insight into the nature of magnetism in these compounds, we 
consider another type of AFM ordering to examine how the relative stability and magnetic 
moments are affected by different AFM ordering. The additional AFM ordering considered 
is a 'checkerboard' type AFM ordering in which the four NN Fe atoms have the opposite spin 
direction to the Fe atom which they surround. This AFM ordering is denoted by AFMl in 
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this paper, and the stripe-type AFM ordering by AFM2. In Table II, the relative energy of 
each AFM type and the magnetic moment on a Fe atom are listed for all the six compounds. 
For each compound, atomic structures optimized in the NM phase are used for all magnetic 
phases to see purely electronic contribution to the total energy differences among magnetic 
phases without structural relaxation effects. 

As shown in Table II, AFMl is more stable than NM phase for all of the compounds, and 
the stability and the Fe local magnetic moment are larger for the antimonides than their 
respective arsenides. Since the AFMl ordering is surely not related to the FS nesting, there 
should be a mechanism other than the simple itinerant magnetism to explain the stability 
of AFMl and its enhancement in antimonides. Furthermore, we find energetic stability of 
AFM2 relative to AFMl phases and magnetic moment in AFM2 phase are enhanced in 
all antimonides compared with respective arsenides, as shown in Table II. This is again in 
contradiction with FS nesting features related to the itinerant magnetism. Therefore, the 
Heisenberg-type magnetic interaction naturally arises as more appropriate description for 
the magnetism in these materials. As the lattice parameters are larger for antimonides than 



their corresponding arsenides, th_e 
reduced band width around Ep 



3d orbitals are more localized as is evident from the 
32l |. Thus the Hund's rule coupling becomes stronger and 
the local magnetic moment is larger for antimonides, as is in Table II. The generally larger 
Fe magnetic moments can explain the enhanced stability of AFMl with respect to the NM 



phase, and AFM2 with respect to AFMl, 
the Heisenberg interaction with J2 > Ji/2 



or antimonides compared with arsenides within 
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29|. 



In the meanwhile, there is clear difference in DOS between AFM2 and other phases calcu- 
lated with the same structural parameters optimized for the NM phase for each compound, 
as displayed in Fig. 3. The NM phase has a finite DOS at Ep, and AFMl magnetic ordering 
does not reduce the DOS at Ep, while it is greatly reduced for the AFM2 ordering. This 
feature indicates that the AFM2 phase involves the ordering-induced FS reconstruction by 
the coupling between the electron and hole surfaces, in contrast to the AFMl phase where 
only the local magnetic interaction is involved. Our result qualitatively agrees with the 



recently suggested model 



3l| in which the itinerant electrons couple to the local magnetic 



moments which are AFM ordered. Even in the case of BaFe2Sb2, where the FS nesting is 
very ineffective as in Figs. 1 and 2, the AFM2 ordering produces the strong perturbing 
potential for the electron and hole bands to be hybridized, resulting in the partial gap in 
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DOS at Ep, as shown in Fig. 3 (d). Other compounds exhibit similar feature in DOS at Ep 
among different magnetic phases, indicating that the presence of partial gap is not sensitive 
to the detailed FS nesting characteristics as it is an induced feature by coupling to more 
robust underlying magnetism of the local moment interaction. 

In summary, we investigate the magnetic properties of known and hypothetical iron 
pnictides by the total energy calculations. We find that our calculated FS nesting feature in 
the NM phase is not consistent with the trend of the magnetic stability that the AFM phases 
are more stable in antimonides than in arsenides. Heisenberg-type local moment interaction 
is more appropriate to explain our results when we consider the larger Fe spin moment found 
in antimonides. Thus our results indicate that experimentally observed stripe-type AFM 
in iron pnictides is mainly driven by local moment interaction, while SDW of the itinerant 
electrons and the partial gap at Ep emerge as an induced order by coupling to the local 
moments. 
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TABLE I: Calculated structure parameters, DOS at Ef (N{Ef)), and Fe magnetic moment (m) 
of LaFeAsO (LFAO), LaFeSbO (LFSO), BaFeaAss (BFA), BaFezSba (BFS), LiFeAs (LFA), and 
LiFeSb (LFA). Both for NM and SDW phases, c lattice parameter is taken as the distance of two 
adjacent Fe layers for easier comparison, zi represents the z coordinate of La, Ba, or Li, and Z2 
represents that of As or Sb. Iron atoms are located at z = 0.5 along the c-axis. 



NM (tetragonal) 

a (A) b (A) c (A) zi Z2 N{Ef) 



LFAO 


3.999 


3.999 


8.706 


0.145 0.640 


1.7 


LFSO 


4.106 


4.106 


9.311 


0.130 0.659 


2.9 


BFA 


3.935 


3.935 


6.314 


0.696 


1.9 


BFS 


4.324 


4.324 


6.315 


0.708 


1.8 


LFA 


3.767 


3.767 


5.967 


0.173 0.734 


2.1 


LFS 


3.995 


3.995 


6.266 


0.211 0.756 


2.6 






SDW (orthorhombic) 






a (A) 


6(A) 


c(A) 


Zl Z2 


m{fj,B) 


LFAO 


5.780 


5.693 


8.875 


0.139 0.654 


2.83 


LFSO 


5.955 


5.844 


9.542 


0.124 0.673 


3.13 


BFA 


5.756 


5.590 


6.520 


0.712 


2.78 


BFS 


6.231 


5.937 


7.246 


0.722 


3.22 


LFA 


5.482 


5.285 


6.190 


0.171 0.745 


2.54 


LFS 


5.830 


5.593 


6.528 


0.199 0.768 


2.95 



FIG. 1: (Color online). Calculated FS of iron pnictides in the NM phase, drawn in the Brillouin 
zone of the \/2 x \/2 supercell (dashed lines) . Hole pockets are represented in blue (dark gray) , and 
electron pockets are in red (gray). In (c) and (d), the conventional simple tetragonal unit cell is 
used rather than the primitive body-centered tetragonal unit cell, for easier comparison with other 
compounds. 
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TABLE II: Stability of magnetic phases and Fe magnetic moments m in fi^ for iron pnictides. 
For each compound, calculations arc done in the optimized structure for the NM phase. Ei is the 
energy of AFMl relative to the NM phases and E2 is the energy of AFM2 relative to the AFMl 
phases, in meV per Fe atom. 



Compound 


El 


m(AFMl) 


E2 


m(AFM2) 


LaFeAsO 


-123 


2.23 


-109 


2.35 


LaFeSbO 


-387 


2.88 


-136 


2.83 


BaFe2As2 


-108 


2.09 


-64 


2.20 


BaFe2Sb2 


-426 


2.80 


-75 


2.78 


LiFeAs 


-45 


1.83 


-99 


1.96 


LiFeAs 


-269 


2.54 


-118 


2.63 



FIG. 2: (Color online). PauU susceptibility Xo(q) for (a) LaFePnO, (b) BaFe2Pn2, and (c) LiFePn, 
normalized by Xo(q = (0)0)) of the arsenide for each type of compounds. 



FIG. 3: (Color online). DOSs of (a) LaFeAsO, (b) LaFeSbO, (c) BaFe2As2, (d) BaFe2Sb2, (e) 
LiFeAs, and (f) LiFeSb calculated in the NM (black), AFMl (red), and AFM2 (blue) phases. 
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